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Abstract 

Recently Kleides et al [IJMPA 11, 1697 (2008)] found a growing rate for magnetic fields 
in ideal plasma cosmologies by making use of general relativistic Friedmann model. This 
growth rate of ^ ~ indicates the presence of a slow dynamo in the universe. More 

recently Hasseein [Phys Plasmas (2009)] have also investigate Beltrami magnetic fields 
in plasma universe. Here general relativistic(GR) MHD dynamo equation, recently given 
by Clarkson and Marklund [Monthly Not Roy Astr Soc (2005)] is used to investigate 
the relation between collapsing of the isotropic universe and dynamo action in ideal and 
dissipative cosmologies. Dynamo action can be supported in these phases as long as 
the kinetic helicity overcomes universe diffusion effects. A cosmological Beltrami flow 
in 3D shows that helicities may act constructively on gravitational collapse and enhance 
dynamo action. A slow dynamo action is found in the static Einstein universe also filled 
with a Beltrami flow. A rotating, shear-free Bianchi type-IX universe, is obtained, by 
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magnetically perturbing the Einstein static model inducing slow dynamos in the model. 

4 

Magnetic field growth of B t, wliicli is stronger tlian Harrison estimate of B ^ fs is 
obtained. CMB limits on the expansion, global rotation and slow dynamos are given and 
a less slower dynamo than the one obtained by Kleides et al, is found with ^ ~ |6|i:. 
PACS: 04.25.Nx,04.40.Nr. key-words:Plasma cosmologies; dynamos. 
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I Introduction 



Importance of the investigation of magnetic dynamos in the GR universe, has earher 
been stressed by Zeldovich, Ruzmaikin and Sokoloff [1]. More recently, Christensson and 
Hindmarsh [2] have considered turbulent dynamos in Early Universe where a homogeneous 
cosmological model was used in connection with string theory. In this model no dynamo 
action was found and a decay of the magnetic was effective. More recently, Barrow 
and Tsagas (BT) [3] in the form of a slow, decaying magnetic field in open Friedmann 
model. In this report one goes a small step further, and computes, based on CM MHD- 
GR dynamo equation, the existence of fast dynamo action [4] in two distinct kinds of 
GR cosmologies. The first is the classical Einstein universe, where it is shown that the 
dynamo action is slow, for a magnetic curvature effect due to magnetic dissipation. Since 
in this case the shear, expansion and vorticity of cosmic fluid vanishes, the CM dynamo 
equation, leads to a simple version even in the presence of non-ideal plasmas. In this 
case, it is shown that, the fast dynamo action can be supported in this cosmology and 
only a slow dynamo is possible. This case represents a cosmological example of Vishik's 
anti-fast dynamo theorem. This problem has been investigated in detail recently by 
Garcia de Andrade [5] in the case of dynamo plasmas [6]. The magnetic field growth 
rate is computed in terms of eigenvalues spectrum of the dynamo operator. In the two 
and three-dimensions, compact Riemannian manifolds, the fast dynamo action has been 
obtained by Chicone et al [7] in diffusive media. This has been done, following previous 
work by Arnold et al [8], in the case of uniform stretching of magnetic field by a steady 
flow in Riemannian space. Besides reproducing the decaying magnetic result of BT in the 
comoving coordinates, fast dynamo action is obtained, when the magnetic field growth 
rate of the Bianchi type-IX is positive. Effective growth rate is defined which includes the 
expansion of the universe. The stretching of universe cosmic plasma is obtained as a source 
for dynamo action. Stretching dynamos by plasma flows, have also been obtained by M 
Nunez [9]. An account of GR cosmological dynamos, is contained in Widrow [10]. Section 
II of the paper, presents the Einstein static universe slow dynamo in diffusive media, 
while in section III, Bianchi type-IX dynamos are presented. In section IV it is shown 
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that a slow dynamo in Bianchi typc-IX model can be induced by magnetically perturbing 
the static Einstein cosmological model. Cosmic global rotation limits on slow dynamos 
is expected from cosmic background radiation (CMBR) data. Section V addresses the 
problem of similar kinetic effects on the existence of dynamo action in expanding phases 
of isotropic universes. 

II Slow dynamos in Einstein static model 

The CM MHD-GR induction equations [11] splitted into {3+l)-spacetimes can be recasted 
in the form used by dynamo theorists 

[1 + e(-e + 2a;)]B = Vx (vxB) + e[AB + 2V(E.a;)- < Ric, B >] + fo (II.l) 
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where A = is the Laplacian, and 

fo = (l + e^e)B-eF- G + Vx(a;xE) - Vx(axE) + ef (II.2) 

F=(a;xB + axE) (II.3) 

where in equation (II.l) 

< Ric,B>= R'jB^ (II.4) 

G= [^ExVe- Vx < (T,E> -Vx(a;xE)] (II.5) 

Here greek letters shall denote space-time indices, while latin letters would represent three- 
dimensional space. Here also, Ric represents the Ricci curvature tensor whose components 
are R^j. Here E, a,e,Q, a and Q are respectively the, the electric held, the acceleration 
in three-space, the constant magnetic diffusion, or resistivity, and the respective Ehlers- 
Sachs optical cosmological parameters, expansion scalar, shear tensor and vorticity vector. 
The later shall be important in the next section when one shall address the rotational 
Goedel-like [12] or Bianchi type-IX universe [13]. one shall not reproduce here the long 
formula for r(a, ©, a, Co, E) and in Einstein static universe in diffusive media it vanishes, 
since acceleration and Ehlers-Sachs optical parameters vanish and the electric field is here 
assumed to vanish as well, which is a reasonable assumption for most dynamo models [1]. 
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Since in the Einstein universe the Ricci tensor is given by the constant Riemann curvature 
contraction, called also Einstein manifold [14] 

Rij = ^9ij (11.6) 

which Einstein latter recognised, as the cosmological constant, one may simplify the above 
CM MHD-GR dynamo equation as 

[1 + €(-e + 2uj)]B = Vx(vxB) + e[AB- < Ric,B>] (II.7) 
3 

This equation can be further simplified, by assuming that 

AB = -curl{curlB) (11.8) 

and that the cosmic flow is a Beltrami flow, where 

curlB ^VxB^^B (II.9) 

which upon substitution into the above equation yields 

AB = -7^5 (11.10) 

Substitution of this last expression and the Ricci equation in Einstein manifold yields 

XB=-e[y^-A]B (11.11) 

since the static case leads to the vanishing of three-speed v. Here the dot represents 
partial time derivative, where one assumes that 

B = Boe^* (11.12) 

Thus the eigenvalue spectrum of dynamo in the Einstein static universe , obtained from 
expression (III. 19) leads to 

X = -£[72 _ A] (11.13) 
Since the slow dynamo is given mathematically by the following limit 

lim,^oRe{\) = (11.14) 



where Re denotes the real part of the growth rate A which is in general complex, though 
here A e ii, where R if the field of real numbers. This is fulfilled by the Einstein universe, 
regardless whether or not, the cosmological constant is bound by the Beltrami constant 7. 
This is due to the fact that the above limit predominates over the fact if the limit above 
is positive or negative. Of course the fast dynamo shall be consider by the positivity of 
the growth rate A. In the next section one shall discuss the fast dynamos in Bianchi type 
IX rotating cosmological model. 
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Ill Bianchi type-IX cosmic fast dynamos in highly 
conducting plasmas 

In this section, one shall consider the Bianchi type-IX cosmology, given by the pseudo- 
Riemannian metric 

ds" = -de + Qi^xY (111.15) 

where the differential basis one-form components are given by 

— —sinx^dx^ + sinx^cosx^dx^ (III. 16) 

= cosx^dx^ + sinx^ sinx^ dx^ (III. 17) 

and 

X^ = cosx^dx^ + dx^ (III. 18) 

In terms of the orthogonal tetrad frame, cu^''^ the pseudo-Riemannian metric maybe written 
as 

ds" = -de + (a;(i))2 + (a;(2))2 + [^i^)f (m.ig) 

wJiere 

a;« = UjX^^^ (111.20) 
In terms of this "dreibein", Uj the pseudo-Riemannian metric can be written as 

9ij = k%j (III.21) 

The Bianchi type-IX cosmological models may expand or contract and have rotation 
according to the Ehlers-Sachs parameters, whose the only one that is missing here is the 
shear. Note that with the help of geodesic equation it is easy to obtain the values for © 

and rotation tensor of components uJik as 

e = - vohk (III.22) 



and 



1 VkVid'^ii 

(jjoi = -iVio = (III.23) 

2 vo 



cuij = -cuji = -^Vkd'^ij (III.24) 
where is the four-spacetime velocity of the cosmic How and 

d^jk — hl^lmnbmj ^Kk ^ (III. 25) 

where 

hj = hij) = + ^iO (III.26) 

and 

hj = hihif^ (111.27) 

By considering extremely highly conducting Bianchi type-IX universe, where the inverse 
of the magnetic Reynolds number Rm, e = RmT^ vanishes, one simplifies the CM MHD 
dynamo equation leads to an effective magnetic field growth rate as 

Ae// = [A + ^e + r-,o] (111.28) 

where one has considered the comoving cosmic frame where — 1 and v vanishes. This 
case is actually similar to the Friedman expanding universe. Since there is in this case a 
hypersurface orthogonality, vorticity vanishes, and only expansion survives. Here F"^^ is 
the Riemann-Christoffel symbol or Levi-Civita connection given by 

r"/37 = la'^'i-degp^ + d^ape + d^go^e] (III.29) 

Computing the above component and the CM MHD GR dynamo equation 

XeffB^ujxB (III.30) 

Since the RHS of this equation, vanishes due to the vanishing of vorticity, one obtains 
that vanishing of X^ff leads to 

A = -[^e + r\o] (111.31) 



which yields finally the equation 



A = i^kk - r"ao] (111.32) 
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for the geodesic triad l^p. The existence of fast chaotic dynamos, depends upon the sign of 
the Riemann-Christoffel component. Thus, let us now compute this component in terms 
of the Bianchi type IX metric as 

r^ao = hk (III.33) 

Now this expression tell us that the growth rale can be given by 

X = -^e (III.34) 

To obtain this result one considers that the magnetic and velocity flow are solenoidal 
divergence-free vectors as 

V.B = (III.35) 

and 

V.v=0 (III.36) 

This result agrees with BT and CM results that, when universe undergoes contraction 
(Q < 0) or stretching in the language of dynamo theorists, there is a stretching in the 
universe; the fast dynamo action is possible since A is positive. In what follows one shall 
prove that this result is possible to be obtained in the less trivial case of the rotating 
Bianchi type-IX, which is the main result of this report. Let us now consider then, the 
equation of the GR MHD-dynamo, in the case vorticity and expansion does not vanish is 

[A + -e + (a.v) + r'^^pv'^jB = Vx (vxB) + a;xB (III.37) 
3 

Note that since the universe is assumed highly conducting e vanishes and the curvature 
effects on the dynamo does not appear in the equation through the term < Ric, B > 

but actually through the Riemann-Christoffel symbol r"^^. This equation can be further 
simplified if one assumes that acceleration a comes from the magnetic Lorentz force as 

a=-vxB (III.38) 

Due to simple vector analysis equations, this yields the following simplifications 

a.v= -v.[vxB] = (III.39) 
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and 

a.B= -[vxB].B = (III.40) 

which reduces the dynamo equation (III. 37) to 

Vx(vxB) +a;xB= [A + ^6 + (a. v) + r"«/3t'^]B (111.41) 

This last equation, can be recasted in dynamo operator form as 

CB = Ae//B (III.42) 

where the dynamo operator can be written as 

£^Vxvx+ujx (III.43) 

while the effective eigenvalue growth rate is obtained as 

Ae//=[A+^e + r^^^;^] (III.44) 

Note that when one is on a syncronised comoving coordinate system where = S'^q, 
expression (III. 44) reduces to the above equation for the Fricdmann dynamo operator. 
Now let us compute the trace of the Riemann-Christoffel symbol, F^q/j contracted with 
the four-velocity as 

T'^apv'^ = Vav" - dav" = e - V. V = e (III.45) 

Since locally the flow is incompressible, or divergence- free. Substitution of expression 
(III.45) into (III.44) yields 

Kff = [A + ^e] (III.46) 

Thus in the case of zero-effective eigenvalue \ff, as a minimum effective eigenvalue, yields 
the following growth rate of the magnetic field as 

A=-^e (III.47) 

which shows that the fast dynamo action may exist in the Bianchi type-IX universe if the 
expansion © is negative. This, to a certain extent, explains why stellar colapsing objects 
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which undergoes this gravitational contraction, may support fast dynamo action. The 
magnetic Geld in this case is 

B' = B^oe~t®* (III.48) 

which shows that when the Bianchi type-IX universe expands, (Q > 0) the magnetic 
Geld decay as in the Friedman case, the role of vorticity being so weak in cosmic dynamo 
here as is in the universe where expansion is much stronger than rotation as measured by 
COBE. 

IV Slow dynamo in Bianchi type-IX from magnetic 
perturbation of Einstein model 

In this section one stablishes a bridge between the last two sections by showing tha 
magnetic perturbation of the static Einstein model from a constant seed magnetic Geld, 
in the absence of shear, may lead to a slow dynamo in Bianchi type-IX model. Let us 
start by the CM equation in the case of highly conductive universe, which is a good 
approximation of our real world. The perturbation scheme is 

B^Bo + Bi (IV.49) 

V = vo + vi (IV.50) 

e = Go + Gi (IV.51) 

uj^ujo + uji (IV.52) 

a = ao + ai (IV.53) 

Since the Einstein universe is static the following simpliGcations arise: All unperturbed 
quantities, but the magnetic seed Geld Bq, vanishes. Substitution of these results into the 
CM dynamo GR-MHD induction equation, yields 

Bi = Vx(vixBo) -a;ixBo = (vi.V)Bo - (Bq.V) vi - cDixBo (IV.54) 

Note that this induction equation is linear in the seed Geld Bq, which allows us to simply 
integrate in time this equation, by assuming that the expansion is almost constant, or in 
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tie small acceleration phase. With the help of the equation 



(IV.55) 



The stretching term is given by 

Bo'ViVi'' = Bo'.iuu" + OiSi''] (IV.56) 

The other term is given by 

[QixBo] = e^jkUJi^Bo'' (IV.57) 

These two terms cancell in equation results into the CM dynamo GR-MHD induction 
equation, yields 

Bi = Vx(vixBo) -tDixBo (IV.58) 



one obtains the equation 



Bi:= SB= -[^QBo + uJixBo] (IV.59) 



which upon integration and Schwarz lemma reads 

l^l<(|le| + l^xb|)t (IV.60) 

where b — ^ is the unit direction of seed magnetic field. Thus since the COBE data, 
given by 

1^1 < 10-' (IV.61) 

provides a bound to the slow dynamo given by expression (IV.60). From dynamo theory 
this is expected, since due to dynamo quenching the dynamo magnetic field cannot grow 
without bounds. So in a certain sense one should say that the COBE data places limits 
not only on expansion and vorticity of the Bianchi type-IX model but also in the dynamo 
action. That is the reason, why the dynamo action is so rare in cosmological large scale 
dynamos, as pointed out previously by Barrow and Tsagas. Other global rotation limits 
from COBE in shear-free Bianchi type model has been proposed previously by Obukhov 
[15] and Garcia de Andrade [16] in Einstein-Cartan gravity, however, from pratical pur- 
poses, rotation in real universes, can be neglected in most situations. Nevertheless from 
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tie dynamo point of view rotation or twist, is fundamental as stretching, and is interest- 
ing to use models that contains rotation to check their influence on dynamo action. Here 
Bi Ki t. By using Zeldovich-Novikov gravitational perturbations in the expanding uni- 
verse, Harrison [17]. has actually found slower magnetic field growth as B ki t^. Taking 
Bi = SB, expression (fV.57) becomes |^| = \Q\t. This result is in agreement with the 
Friedman model result |^| = |^|t ~ 10~^. 

V Dynamo action in expanding universes 

Based on their MHD dynamo equation, Clarkson and Marklund [11] have shown that 
the cosmological gravitational collapse of the universe shall certainly leads to dynamo 
action while in expanding universes, the magnetic field would decay as shown in BT 
paper refereed above. In this section one shows that the universe expansion phases can 
actually support dynamo action as long as the mean electromotive force, given by 

E^vxB^aB (V.62) 

and a kinetic helicity effects overcomes universe expansion 0. This simple proof can be 
done by substitution of expression (V.62) into the dynamo GR MHD formula (H.l) by 
neglecting electric fields, yields 

B = Vx (aB) + e([AB]- < Ric, B > +fo) (V.63) 

In the case of Beltrami flows, equations (II.9) and (11.10) yields 

B(l + ^GA) = [{a - 6)7 - (^e + h{E -3p + AA+ '^e'))]B (V.64) 

By assuming the diffusion is small and 0{e^) — > 0, one is able to obtain the constraints 
to helicity spectrum in order to obtain the dynamo action without necessarily to have 
collapse or G < 0. Here the relation 

i?,, = ^A(£; + A-^e2) (V.65) 

between the Ricci tensor, and energy density e, and the cosmological constant A is used. 
Here hij is the projective 3D metric. By assuming that a marginal dynamo is not a 
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dynamo action, one may find a simple degenerate spectrum of expansion by assuming 
tliat B vanislies on tlie LHS of induction equation (V.65) as 

[(a - e)7 - (^e + ^X{E - 3p + 4A + ^6^))] = (V.66) 

Tliis second order algebraic equation in the universe expansion yields the real helicity 
spectrum as 

o? - 2ea + ^ (£; - 3p + 4A + ^e^) = (V.67) 
which yields the following constraints for real spectrum 

{R + 2A+^e^) <3p (V.68) 

which yields a highly pressured and early phase of the universe where helicity acts. Note 
that since the expansion appears squared in this condition does not mind if the universe 
is expanding or contracting for marginal dynamos. Note, however in this case, unless the 
Ricci scalar and cosmological constant be negative relation (V.68) is violated. But this 
is exactly the gravitational collapse case. However, in the infinite conductivity 
simple proof of the fact that even in expanding phases of the universe helicities could be 
strong enough to enhance dynamo action can be given. This can be done by considering 
that the above equations in this section reduce to 

X = aj--e (V.69) 
3 

This shows that a fast dynamo action is possible even when the expansion of the universe 
(6 > 0) is taking place as long as the product of the kinetic and magnetic helicities are 
positive and strong enough to overcome expansion. This happens in the early phases 
in the universe. Of course this shall be stronger enhanced when the contraction phases 
taking place. That could explain why BT results where even in slow expansion a decay 
of magnetic fields is possible. Thus the universe expansion phases are more willing to 
possess decay magnetic fields than dynamo action. 
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VI Conclusions 

Thus one may conclude that by the analogous stretch-twist-fold Vainshtein-Zeldovich 
mechanism [18] in cosmology given by expansion-rotation-curvature, one is able to ob- 
taining amplification of the magnetic field, from a seed magnetic field by dynamo mecha- 
nism, producing at least a slow dynamo. Therefore the slow dynamos obtained here, from 
magnetically unstable Einstein static universe, endowed with a seed field, seems to be 
in agreement with previous ideas that the seed fields generated during phase transitions 
in the radiation era, have typical coherence lengths that, by its smallness,may destablize 
the dynamo. As also pointed out by Barrow et al [19], fields which survives inflationary 
epoch, are in general too weak to sustain a fast dynamo. It is interesting to note that, the 
smallness of the global rotation of the universe seems to lead to a slow dynamo, which 
eventually decays upon expansion. This eventually does not happen in the case one intro- 
duces mean held helicity into Clarkson Marklun MHD dynamo equations, as long helicity 
overcomes expansion. 
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